function [x,y]=betapdfcheck(para1,para2);

close all;

a = (1-para1)*para1^2/para2^2 - para1;
b = a*(1/para1 - 1);

x = 0:.01:1;
y = betapdf(x,a,b);


fmt = '%0.5g';
plot(x,y);
title(['Beta with (' num2str(a,fmt) ', ' num2str(b,fmt) ')']);
xlabel(['\mu = ' num2str(para1,fmt)  '    \sigma = ' num2str(para2,fmt) ]);
